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Abstract 



We demonstrate, to next-to-leading order accuracy, the cancellation of the infrared 
singularities in the impact factors which arise in the QCD description of high energy 
processes A + B ^ A' + B' of colourless particles. We study the example where A is a 
virtual photon in detail, but show that the result is true in general. 



1. Introduction 

There continues to be much activity in the study of semihard QCD processes, that is pro- 
cesses at small values of x = Q'^ j s where is a typical virtuality and ^/s is the centre-of-mass 
energy. Many intriguing low x phenomena have been observed (see, for example, ref. since 
the announcement of the sharp rise seen in the proton structure function at HERA at x de- 
creases. The theoretical understanding of these phenomena is a non-trivial task. One important 
type of process is high energy diffractive qq (or vector meson) electroproduction, which to lowest 
order proceeds via two gluon exchange. Much data are becoming available for such reactions. 
However the size of the QCD radiative corrections is not known. There are indications both 
experimentally and phenomenologically that these corrections are very important 0. Another 
relevant process is the total cross section for 7*7* collisions. 

The most common basis for the description of small x processes is the BFKL equation 0. 
In this approach the amplitude for the high energy process A + B A' + B' a.t fixed momentum 
transfer may be symbolically written as the convolution 

$A'A ® G ® ^B'B, (1) 

see Fig. 1 and eq.(^ below. The impact factors $a'a and $b'_b describe the transitions A —>■ A' 
and B ^ B' shown by the upper and lower blobs of Fig. 1, while G is the Green's function for 
the two interacting Reggized gluons. This representation is valid both in the leading logarithmic 
approximation (LLA) , when only the leading terms {as In s)" are resummed, and in the 
next-to-leading approximation (NLA) 0, when the ^^(as' In s)" terms are also resummed, 
not only for forward scattering Q with t = 0, but for the non-forward case as well. The 
explicit form of the kernel for forward scattering (with colourless exchange in the t channel) 
was found in ref. |]. A number of subsequent papers (see, for instance, [0]) were devoted to 
investigations of its consequences. 

Both the kernel of the integral equation satisfied by the Green's function and the impact 
factors are unambiguously defined @, |, ||] in terms of the gluon Regge trajectory and the 
effective vertices for the Reggeon-particle interactions, which have been calculated in the NLA 
in a series of papers P|-|T^. The cancellation of the infrared singularities in the kernel was 
explicitly demonstrated in ref. [Q. We stress that, for brevity, we are using the terminology 
"infrared singularities" to include both soft and coUinear singularities. Strictly speaking the 
Green's function for the Reggeized gluons, taken by itself, is not physical, since only the am- 
plitudes describing the interaction of colourless particles are measurable. It is these physical 
amplitudes that must be free of infrared singularities. However since it has been shown that the 
kernel of the equation for the Green's function is infrared safe, it should therefore be possible to 
construct impact factors which are infrared safe. Recall that we are discussing impact factors 
for transitions between colourless particles and, consequently, about a colour singlet state in 
the t channel. 

^The impact factors were defined in refs.[|[ ^ Q in a different way; see tlie comment in Section 7. 
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Here we study such impact factors at NLA accuracy. We demonstrate explicitly the can- 
cellation of the infrared singularities in the impact factor of the virtual photon and, moreover, 
give a proof that the impact factors describing the transitions between colourless particles are 
infrared safe in general. We concentrate on the virtual photon because of the importance of 
the NLA for diffractive qq electroproduction, as mentioned above. Our study may therefore be 
considered as a first step in analyzing the next-to-leading radiative corrections to this experi- 
mentally accessible process. 

The outline of the paper is as follows. In Section 2 we give a detailed definition of the impact 
factor, concentrating on the example of that for a virtual photon. In Section 3 we calculate 
the singular contribution which comes from the virtual correction to the LLA approximation 
for the effective virtual photon-Reggeon interaction. Then in Section 4 we consider the effects 
arising from an additional real gluon emitted in the 7*-Reggeon interaction. In Section 5 we 
draw these results together and demonstrate the cancellation of the infrared singularities in the 
impact factor of the virtual photon in the NLA. A proof of the cancellation in the general case 
is given in Section 6. Finally, in Section 7 we present our conclusions. 



2. The definition of the impact factors 

First we make explicit the convolution (|I|) for the amplitude T describing the process A -\- 
B ^ A' + B' shown in Fig. 1. To be specific we calculate the s-channel imaginary part of the 
amplitude 



q\q{ J qi^qi 

(2) 

duj \f sy - 
— '^uj[qi,q2'^q) 

Sn/ J 



(27r)^-2 Js-ioo 2m 



where the momenta are defined in Fig. 1. For convenience we introduce q[ = qi — q, where 
g ~ is the momentum transfer in the process A + B ^ A' + B' 

q = Pa-Pa' = Pb' - Pb- (3) 

^A'A and ^B'B are the impact factors which describe the A A' and B ^ B' particle 
transitions in the particle- Reggeon scattering processes shown in Fig. 1, and G^^ is the Mellin 
transform of the Green's function for Reggeon- Reggeon scattering. We emphasize that the zig- 
zag intermediate lines in Fig. 1 denote Reggeons, and not gluons — the Reggeons would only 
become gluons in the absence of interaction. The integrations in (H) are performed over D — 2 
dimensional vectors which are transverse to the momenta pa and Pb oi the initial particles. 
The space-time dimension, D, is taken to be D = 4 -|- 26 in order to regularize the infrared 
divergences. 

Representation (H) is derived assuming that the energy scale sq is independent of q^ and 
5 or at least has a factorizable dependence on q^ and Then any sq dependence induced 
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by using sq as the energy scale in can, to NLA accuracy, be transferred to the impact 



factors. This was shown in ref. [|T^] for forward scattering; for the non-forward case the proof 
is essentially unchanged. 

We consider a colour singlet state in the t-channel (or, to be more precise, vacuum or 
Pomeron quantum number exchange), since only in this case will the infrared singularities 
cancel. Thus we can suppress the colour group representation indices of the impact factors and 
the Green's function. We also note that representation applies not only to elastic processes, 
but also to reactions where A' and B' may each correspond to a group of particles with invariant 
mass independent of s. Finally the normalisation adopted in (0) is that of which seems to 
be convenient for non-forward processes. It differs at t = from that used previously such that 

GMi,q2\q)\q=o = <llqlGio{qi,q2) (4) 

where G^{qi, is the Mellin transform that was used for the forward case [Q. 

We may express both the kernel of the equation for the Green's function and the impact 
factors in terms of the gluon Regge trajectory and the effective vertices for the Reggeon-particle 
interaction . We begin with the high energy amplitude T^^^ for the process A + B ^ a + b 
of Fig. 2, which is an intermediate process of our reaction A + B ^ A' + B'. We have colour 
octet Reggeized gluon exchange in the t/j channel (tji = g|, ~ qj^) with negative signature so 
that 

where F^^ and F^^ are the Reggeon vertices for transitions A —>■ a and B —>■ b and 1 + u!{t) is 
the gluon Regge trajectory. Note that this representation is true for particle-parton as well as 
for parton-parton processes. By parton we mean a coloured quark or gluon, whereas a particle 
is colourless. In the parton-parton case the vertices are known in the one-loop approximation 



rO| . One of our aims is to identify the particle-parton vertices in the same approximation. The 
trajectory function, 

j{t) = 1 + uj{t), (6) 

is independent of the nature of the scattering objects, and is now known in the two-loop 
approximation |Tl|. The summation over the repeated index i is over the 8 coloured gluon 



states. 

In the LLA the impact factors are given by 

^%iqn,q) = E / ^ rl^fe) F^^lg;,)* rfp. (7) 

where the number of colours Nc = 3 for QCD and F'^'^^ means that the vertices are evaluated 
in the LLA or Born approximation. The sum {a} is over all states a which contribute to the 



3 



A ^ A' transition, and over all discrete quantum numbers of these states. The phase space 
element dpa of a state a consisting of particles with momenta is 

dPa = (2vr)^ 6(pA-qn-i: li (n^^u-tj, ^ (8) 

\ nea / nea l^'' J ^-'^n 

while the remaining integration in is over the squared invariant mass of the state a 

Ml = (pa - QR? = (Pa - Qr?- 

In the LLA the impact factors do not depend on sq. In fact the energy scale is arbitrary in this 
approximation. 

To be specific, we choose the relevant example of the impact factor of a virtual photon 7*. 
In order to define the vertices we may study the high energy interaction of the 7* with 
any particle. Clearly we should take the simplest choice for particle B — a massless quark 
with momentum ps = P2 and P2 = — and introduce the light-cone momentum pi so that the 
virtual photon momentum is 

Q2 

Pa = Pi —P2 and s = 2pi.p2, (9) 

where —Q^ is the virtuality of the photon. 

The lowest-order 7*g process which has a non-vanishing cross section at large s is 7*g — >■ 
{qq)q where the qq pair is produced in the fragmentation region of the photon. The Feynman 
diagrams are shown in Fig. 3. We may evaluate these diagrams using the trick of writing the 
metric tensor in the gluon propagator as 

^ + + ^M^, (10) 

s s 

and retaining just the first term. Only this term gives a leading (growing as s) contribution 
to the amplitude. In this way we arrive at the 7*g {qq)q Born amplitude in the form of (j^) 
with the lower vertex (describing the q ^ q' transition via the Reggeon field) given by 

r??: = guip'^) ^ f u{p2) = g{q'\f\q)5,,y^, (11) 
and the upper vertex (describing the 7* — > qq transition) given by 

rSl;*(gH) = (g|tig)r^,,^.(gH); F,,,,. (g^) = + F, (12) 

with 

r. = '-^uiU) + vii_) (13) 

F, = - ^ uH^) + (14) 

S ^-i- 



The components and correspond to the upper vertices in diagrams 3a and 3b respectively. 
The notation used in (|1T])-(|1^) is as follows: A denotes the quark helicity, g is the QCD coupling 



and (a|t*|6) are the matrix elements of the colour generator in the fundamental representation. 
Also ecq is the electric charge of the quark (in the qq pair), e^* is the polarization vector of the 
photon, i± are the q and q momenta shown in Fig. 3 and, finally. 



^ + x±Q2. (15) 



We have taken the quark to be massless and have used the decomposition 

i± = x^pi + + 4±, 4 = (16) 

sx± 

In fact the expressions (p!TD -(p!4D for the vertices could have been written immediately by noting 
that the vertex factors^ igf'J)i^2/s describe the interaction of a Reggeon with quarks having a 
large component of momentum along p2,i respectively. In other words in the LLA the Reggeon 
acts as gluon with polarization vector —p2,i/s when it interacts with a quark having a large 
component of momentum along pi 2 respectively. Taking this into account we can represent the 

.(0)i 



vertices F^^^^, by the diagrams of Fig. 4. 

The vertex F^j*^, , (g^) can be calculated in just the same way. It is given by (|T2|)-(p!4D with 



replacements 7* — > 7* , qn ^ qR, ^± — > A'^, where 



and —Q'^ is the virtuality of the final photon. 
Let us denote 

/ N p(0)i 



A'^= ^^^±^ + x^Q'\ (17) 



iV2 -1 



E rS!r(?«) ^'Z^'ilRT = U^±,£±,<1r,q'r)- (18) 

c ^ {gg} 

Besides the variables explicitly shown in (|18|), the function fgg also depends on the fixed pa- 
rameters Q^, Q' ^ and on the photon polarisation vectors. The photon impact factor in the 
NLA is given by the integral of the function fgg over 

^ dpgg = 25(1 -x+- a;_)(27r)^-i5(£+ + + q^) 



27r 



dx+ dx- d^-H+^ d^-H-± 
^ 2^7 2x_ (27r)^-i (27r)^-i ' ^ ^ 



^These are the simplest vertices in the effective action for Reggeon interactions ||T^ which contains all LLA 
Reggeon-parton vertices. 
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In fact this impact factor differs by only trivial factors from the analogous one in QED, which 
was obtained for some particular cases many years ago [0. The explicit form of the function 
fqq is not important for our purposes. The crucial properties of this function are: 

fqq{x±,£±, 0, q'^) = fgg{x±, £±, q^, 0) = 0. (20) 

These properties can be verified directly by using expressions (|l^)-(0) for T^^^!^. (g^) and corre- 
sponding expressions for T^^^.l^„{q'j^). But, in fact, these properties have a quite general nature 
and are consequences of gauge invariance. These properties guarantee that the impact factors 
of colourless objects satisfy the equations 

<I>a'a(0, q) = $A'A(q, q) = 0, (21) 

that are necessary for the corresponding cross sections to be finite. The next important fact is 
that the integral of the function fqq over the measure ([T9| ) is not singular, so that we can take 
x± ~ 1, £-± ^ /i^, with some nonzero /i depending on Q^, Q' ^, q^, {qj^ — q)^. For our purposes 
we can consider all these parameters to be of the same order. In the following we use /z^ to 
denote this order. 

In the NLA the expression (|^) for the impact factor has to be changed in two ways. First 
we have to take into account the radiative corrections to the vertices F considered above. We 
note that the corrections depend on an energy scale. In the definition (^) the energy scale was 
put equal to —tji. Since we calculate impact factors for scale Sq we must allow for the change 
of the vertices arising from the change of scale. It is clear from (|^) that the ratio of vertex 
factors evaluated at scales si and S2 is equal to (■S2/si)'^*-*^''/^ for a Reggeon of momentum qr. 
Secondly, in the sum over {a} in (|^), we have to include more complicated states which appear 
in the next order of perturbative theory, such as states with an additional gluon. For the photon 
impact factor qqg is the only such additional state. But here we have a problem. The integral 
over M^, which is well convergent in the LLA case, becomes divergent when an extra gluon is 
included in the final state. The divergence arises because the gluon may be emitted not only in 
the fragmentation region of particle A, but also in the central region and in the fragmentation 
region of particle B. All three contributions have to be carefully separated and the last two 
assigned to the Green's function G and the impact factor ^b'b respectively. The separation 
was described in ref. IHl. The result for the non- forward case is H 



1 ^ f dM^ 

1 f d^-\n' 



Jn% - 1 ^ 27r 



2 J qU% 



Ki'\qn',qR; q) in 7 (22) 

\{qR'-qRrsQj 



+ qUl' S{qR - qR.d [^{tR) In (^^j + u;(4) In (^^ 

where the intermediate parameter s\ should go to infinity. The second term on the right 
hand side of (p^) is the subtraction of the contribution of the emission of the additional gluon 
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outside the fragmentation region of particle A. The dependence on sa vanishes because of the 
cancellation between the first and second terms. K!^^^ is that part of LLA kernel of the equation 
for the non-forward Green's function which is related to real gluon production M 



(23) 



The kernel describes the transition from the t channel state of two Reggeons with momenta qr 
and — = q — into the state with q^i and — g^/ = q — qRi. The terms which involve the 
Regge trajectory in (|22| ) arise from the change of energy scale in the vertices F. The trajectory 
function uj{t) has to be taken in the one loop approximation 

where F(x) is the Euler function and 

f ^ g^NcVil - 6)/(47r)^/2. (25) 



As stated above, to evaluate the photon impact factor in NLA we have to evaluate the 
virtual corrections to the 7* — > qq vertex of (|l^) and to include the effects of the emission of 
an additional gluon (that is to find the 7* — qqg vertices and to include their contribution in 
(|22|)). We begin with the virtual corrections. 



3. One-loop corrections to the 7*i? qq vertex 

One way to determine the corrections to the 7*i? — > qq vertex is to calculate the one-loop 
corrections to the diagrams of Fig. 3 with a colour-octet state in the channel. We keep not 
only the In s terms, but the constant terms as well. We will find the coefficient of the In s term, 
which is due to gluon Reggeization, is 

2s 



rgtC^Xt^), (26) 



t 



R 



as it must be. If we compare the calculated corrected amplitude with representation (^ and 
use the knownQ expression for the lower vertex F*; ^ to one-loop accuracy |rD|, then we may 
extract F*^,^. in the one-loop approximation. 

In practice we can determine the 7*i? — > qq vertex to one-loop accuracy in a more economical 
way by focusing attention directly on the vertex and not on the amplitude. To be specific we 
work in the Feynman gauge. There are several ways in which we can add the extra virtual gluon 
line to the diagrams of Fig. 3. It is clear that if one end of the extra gluon line is attached to the 
g or g of the upper gg pair and the other end is either again attached to one of the upper pair 
or to the exchanged gluon, then the correction is to the upper vertex. Similar considerations 



^It was determined using the same procedure for the quark-quark scattering process. 
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apply to the lower vertex. It is also clear that half the gluon self-energy correction has to be 
included in the upper vertex and half in the lower vertex. The only uncertainty comes from the 
diagrams of Fig. 5 with two gluon lines exchanged between the upper and the lower vertices. 
Perhaps a better terminology, which we shall adopt, is to say two gluons exchanged between 
the upper and lower blocks. We analyse these two-gluon exchange diagrams first. 

3.1 Two-gluon exchange diagrams 

We start with diagram 5(a). The basic integral for this diagram is 

J ^ f I (07) 

It is convenient to use the Sudakov decomposition for the gluon momentum, 

k = Ppi + ap2 + k±, d^k = ^dadpdk^'^, (28) 

and to consider three regions of the variables a, /3 : 

the central region |a| < ao, \&\ < Po 

the A region > /3o (29) 

the B region \a\ > ag, 

with CKo and Po chosen so that 

ao < 1, /3o < 1 and sao/^o > A*^- (30) 

Recall that we use n to denote the order of magnitude of typical transverse momenta. In fact 
we will take the limit saof^o ^ cxDj while ao ~^ and f3o 0. We see that regions A and B 
overlap. The whole a, j3 plane is covered by the sum of the three regions, with the subtraction 
of the overlapping subregions where both of the variables are "large", that is |a| > ao and 
1/9 1 > /So- If we decompose the momentum transfer 

(IR = (^rPi + OiRP2 + qR^ (31) 

then we have ^ ^ 

sI3r ~ q% saR -Q^ - ^ - ^. (32) 

Combining this with the fact that all transverse momenta are of order n, it is easy to show 
that the overlapping regions where |a| and \f3\ are both "large" give a contribution of relative 
order < fi^/saoPo, which vanishes in the limit s oo. Thus I^a is given by the sum of the 
contributions from the three regions. The result for the central region {\a\ < ao, \(3\ < (3o) is 



.central ^ 2r(l - s) T^s) {qlY 

5a (47r)^/2 r(2e) 2sx_ 



In (^-^^^ + - i'il - e) + 2V(e) - 2i,{2e) 



(33) 
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where ip{x) = T'{x) /T{x). We emphasize the important fact that the contribution of the central 
region is universal (i.e. process independent). The above equation gives not only the singular 
1/e term, but also the general e dependence. For completeness we give the contributions from 
regions A and B which are non- vanishing as e — 0: 



-'5a 



2T(1 



(47r)^/2 2sx. 




where 



and 



L I 



fif 

L{x)= / -ln(l-t), 
Jo t 



In 
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In' 



(34) 



'5a 



2T(1 



{<!%) 



(47r)^/2 




+ ^(2e) 



(35) 



r(2e) 2sx- 

We checked by an independent calculation of the integral of (0), using Feynman parameters, 
that is indeed given by the sum of the three contributions, and hence that there is a negligible 
contribution from the regions where both |a| and are "large". 

We now study how the basic integral ( p7|) enters the matrix element. We use the same trick 
that we used before for the metric tensors in the gluon propagators, see (plOl) . In this way we 
find factors 

M-^--m2 and MA- mi (36) 

respectively in the numerator for the upper and lower blocks in Fig. 5a. It is clear that we can 
neglect ^ in the central region. Since we can omit the contribution from the region where both 
|a| and \P\ are large, we are also able to neglect ft in the numerator of the upper (lower) block 
in the A{B) region. We can then factorize from the numerator the corresponding term which 
arises in the amplitude of Fig. 3a. To see this we use the anticommutation relations of the 7 
matrices to move l^-ij)'^) to act on the spinor f (£_)('u(p2)), and then use the Dirac equation. We 
also have factorization in colour space. We require negative-signature colour-octet exchange 
in the t/j channel, and so we need the antisymmetric colour-octet projection of the two-gluon 
state. The colour structure is thus identical to that of the diagrams of Fig. 3. 

Therefore it turns out that the contribution of the "central" region of Fig. 5a to the 7*g — >■ 
(gg)g amplitude T*^^) of (|^) is 



central 
5a 



f (g|f|g) r.r(?)Xt«) 

^R 



— s 



In 1^-^ J + 0(«o) + 0(/3o) 



(37) 



where is the 7* — > qq Reggeon vertex of (p!3|), FL is the q ^ q' Reggeon vertex of (|n 



is the one-loop contribution to the trajectory function of (|25|) , and the (f) functions are 

1 



Inz 



[(^(1) - ^(1-e)] - ^(e) + ^(2e). 



(38) 
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Diagrams 5b, 5a' and 5b' can be computed in exactly the same way. The contribution of the 
central region of Fig. 5b is given by (^) with the replacement of Tf, of ( ]T^ ) for F^; while the 
contributions of 5a' and 5b' are obtained from those of diagrams 5a and 5b by the replacement 
s —s and by changing the overall sign (due to antisymmetry in colour space). Collecting all 
these four central region contributions together we therefore have 



T 



central 



.(0)i 



R) 




0(ao) + 0(/5o) 



(39) 



where we have used (|T2p. We see that the logarithmic part of the contribution coincides with 
the first term in the expansion of (|^) in uj{tR); that is, it is responsible for gluon Reggeization. 
The remaining pieces, 0(ao) and ^(/Jo), must be included as corrections to the corresponding 
vertices. The important fact is that contribution ( |39D of the central region is universal. That 
is this result for 7*g — > {qq)q applies to any process A + a ^ i? + 6, as in Fig. 2, with the same 
expression in square brackets, provided that we use the corresponding vertices. 

We now turn to the A region (|/3| > (3q, see (^)). We first consider the combined contribu- 
tion of the diagrams of Figs. 5a, 5a' and 6a. To evaluate diagram 6a we again use the familiar 
trick ([T0|) for the metric tensor in the propagator of the gluon with momentum qji- In this way 
we show that we can factor out the vertex Fg?^. The same factorization applies to diagrams 5a 
and 5a' in the A region. Moreover here we can replace 

(p'2 - kf + = -s{l3 - /?k)(1 - a - an) - {q^ - kf + ^ -s/3 
{P2 + kf + ^0 = s/3 - + iO ^ s(3. 



(40) 



Then taking into account the difference in signs for diagrams 5a and 5a' due to colour anti- 
symmetry, we see that they give equal contributions in the A region. After the factorization 
of the common factor F^?^ we find that the remainder does not depend on the properties of 
the lower block (vertex). Thus the sum of the remaining contributions of diagrams 5a, 5a' and 
6a must give corrections to the vertex Fg^^^*. Moreover this sum is given by Fig. 7a, where we 
have introduced the off-mass-shell gluon-gluon-Reggeon vertex of Fig. 8 



liabikl,k2) 



ab 



-g'"'P2-{k2-k,) - pi^{2k, + k2) 



(41) 



+ p^(2A;2 + kif - 2{ki + k2 



V2P2 



P2 ■ {ki - k2 



where are the generators in the adjoint representation in colour space; 

(0)i 



cab' 



In the same way we factor off the lower vertex F^, ^ from the A region contributions of 



diagrams 5b, 5b' and 6b, and find that their sum is given by Fig. 7b. It is also clear that the B 
region contributions of diagrams 5a, b, a', b', in which the upper vertex F^^^^. has been factored 
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off, will give the corrections to the lower quark-quark-Reggeon vertex T^/ g. We are, however, 
not concerned with these latter corrections here. 

Now we turn to the remaining two gluon exchange diagrams, 5c, c'. Unlike the previous 
diagrams here both the produced q and q interact, and they contribute in quite different regions 
of the Sudakov parameters a and (3. Only the region of small a, \a\ ~ fs^/s, is important. The 
reason is that the two virtual fermions in the upper block have momenta components of order 
unity along pi, and ~ a + 0{fi'^/s) along p2. The denominators of their propagators are 



SX- 



X ^{x^ + (5r -(5)\aR-a^ ^ j - (£+ -V q^- k)' + tO j , (42) 

whereas we find the numerators (after performing the trick with the metric tensors in the gluon 
propagators) do not contain any a dependence. Therefore the integral over a converges with the 
main contribution coming from the region |a| ~ yU^/s. Let us now consider the f3 integration. 
The contribution of both diagrams 5c and 5c' to the colour-octet, negative- signature amplitude 
in the channel contains a factor, from the fermion propagators of the lower block, of 

' ' . J . (43) 



V/3-fcVs + iO (3- (3R + {qR-kf/s-iQ^ 

In the limit s — oo it gives 2P(1//?), since (3r ~ q\/s, where P means the principal value. 
For |a| ~ yU^/s, the (3 dependence in all the other propagators can be neglected at small (5. 
Thus the region \(3\ < (3q does not contribute. Hence diagrams 5c, c' contribute only in the 
A region. As before, we can combine them with diagram 6c and factorize off the vertex vf, 



(0)i 



? .9 



The remainder can be represented by diagram 6c with the off-mass-shell gluon-gluon-Reggeon 
vertex of (^. 

Let us summarize the analysis of the two-gluon exchange diagrams of Fig. 5. These diagrams, 
which are responsible for gluon Reggeization, also contribute to the Reggeon vertices. The 
contributions to the 7* —>■ qq vertex, which we are interested in, come both from the central 
region {\a\ < ao, < /3o) and from the A region (|/5| > j3o). The contributions coming from 
the A region can be combined with the contributions of Fig. 6 and represented by the diagrams 
of Fig. 7 with the gluon-g luon- Reggeon vertex given by (|41|). The contribution coming from 
the "central" region can be extracted from (^) 



^^.(central) ^ ^.{qlT ^^a* 0(/?o), (44) 



where we have used (0), and where (f){z) is given by (pH]). The intermediate parameter (3q 
cancels when we combine (|4^) with the contributions of Figs. 7a,b. (The contribution of 
diagram 7c does not depend on Pq.) 
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In addition to the two-gluon exchange contributions, we have also to include corrections to 
the 7* qq vertex coming from the diagrams of Fig. 9, as well as from the quark self-energy 
diagrams, and also from half of the self-energy contribution of the gluon of momentum g^j. The 
latter contribution is not infrared divergent and will not be discussed further. 

3.2 Loop diagrams within the produced qq pair 

The infrared divergent contribution of the quark self-energy diagrams can be found as 
follows. First the self-energy of the internal fermion lines is not infrared divergent. For the 
external lines, in dimensional regularization, the ultraviolet and infrared divergences cancel. 
Therefore it is simplest to determine the infrared divergence by taking the ultraviolet divergence 
with the opposite sign. It gives 

A -piCquark self-energy) _ -p{0)i ( 9 dp \ s 

- ^ [--J^)- ^^^^ 

It remains to calculate the divergent contributions of Figs. 7 and 9. Unfortunately, in 
the Feynman gauge, there are many artificial divergences in the contributions of the separate 
diagrams, which cancel in the sum. Although calculations can be done more easily by the 
choice of an appropriate gauge, we prefer to work in the Feynman gauge in order to escape 
ambiguities of the light-cone gauges. 

Our procedure is to calculate Z)-dimensional integrals over the virtual gluon momentum 
k using the Sudakov components. We first perform the integration over a using the residue 
theorem. We then introduce an intermediate parameter xq to consider the coUinear singularities 
in the regions < Xq and > Xq separately. We take Xq small, Xq -C 1. We let it tend to 
zero but not before e ^ 0, so that eln(l/a;o) ^ as e ^ 0. On the other hand we take jSo —>■ 
before e — 0. Recall that all calculations at next-to- leading order are performed by first taking 
s — oo and then e — 0. 

For \P\ < Xq all three diagrams of Fig. 7 contribute, as well as diagrams al and bl of Fig. 9. 
At first sight it may seem strange that the diagrams of Fig. 7 are divergent in the 'soft' region. 
However recall that they contain the gluon-gluon-Reggeon vertex and receive contributions to 
the 7* — >■ qq vertex from the diagrams of Fig. 5. The contribution of diagram 7a is 



A -p (soft 7a) _ _ 2 



qlr In ----- In 



r(2e) /5o e2 e \xqx_ 



Ta, (46) 



while that of 7b is obtained by the substitution q <->■ q, that is by F^, A_ ^ A+ and 

X- <-> x+. The ln(l//3o) terms cancel the exactly analogous terms coming from the central 
region, see (|^) and (|38D , as they must. Diagram 7c, which is symmetric with respect to the 
q q substitution, is given by 



A p (soft 7c) ^ _ 2 




Ta + iq^q). (47) 



12 



Turning now to diagram 9al we have 



Ar 



(soft 9al) 

mn* 



2 1 

9 ^ 



N 



c 



^ + - In 







X+X- 



(48) 



while the contribution of 9bl is obtained by the substitution q q. The total contribution 



from the central region [ (|4^) , (pq) ] and the soft region [(|46| ) - (|48|) ], plus the q q diagrams of 
7b and 9bl] is 



AT 



i(soft + central) 



N6 



.(0)i 



(49) 



We now study the diagrams with coUinear singularities. As mentioned above, in the Feyn- 
man gauge there are artificial singularities in the separate diagrams, which cancel each other in 
the sum. The final answer is very simple and could be written down straightaway. Nevertheless 
we check explicitly that the cancellation indeed occurs. 

For instance diagram 7a has a singular contribution from the region in which k is coUinear 
with i_ 



A pi(7a coUinear) 
9917* 



with the non-factorizable contribution 



In 



X- 



Xo 



1 (g|f|g)r, - Rl 



'1 — z)x- 



e Jo 2[qR + zi^)^ 



(50) 



(51) 



It is easy to understand the origin of such non-factorizable terms, and to follow their mutual 
cancellation using the gauge properties of the gluon vertices. The gluon-gluon-Reggeon vertex 
(|1|) gives 

(52) 



^g rpi 

s 



kfp^ + {k2-P2)K 



since p2 ■ ^2 — —P2 ■ ^1 as the Reggeon momentum qn = ki — k2 is almost transverse; thus p2 ■ qn 
does not grow with s, whereas p2 ■ ki and p2 ■ ^2 are proportional to s. A vertex describing 
the interaction of a coUinear gluon and a fermion gives us a gluon momentum which has to be 
contracted (via the gluon propagator) with another gluon vertex. For example, in the region 
where k and £_ are coUinear, k ~ — the matrix element of diagram 7a contains the factor 



(/_ + l^)Yv{L 



2(1 - z)r t;(£_) 



-2^^ — ^ri;(£_). (53) 



Now we use the gauge properties of the vertices. For diagram 7a we put k2 = k and ki = q^ — k 
into the gluon-gluon-Reggeon vertex of (|5^). The first term gives a contribution proportional 
to Ta, and the second term the non-factorizable contribution i?^ of (|5lD, where we recognize 
the gluon propagator with momentum qn — k with k = —zi-. 
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The contribution of diagram 7b has a form analogous to (|50[) , coming from the region where 
k is coUinear with 

^^.(rbcouinear) = _ ^ (^1^ ^ _ (^q\f\q)T, - Rl (54) 

where 

Rl ^ -e^^ /' (55) 
e Jo 2(g^ + 2;^+)^ 

Diagram 7c has contributions both from the region where k is colhnear with and from 
the region where qn — k is colhnear with i_. It has the non-factorizable form 

^^^Jccollinear) ^ ^ ^ ^ (gg) 

where R\ and 

= M±±iiL±it|r„(,_, (57) 

e 7xo/x_ 2; ^ +^ s(£+ + gjj^ + z£_)2 ^ ^ ^ ^ 

come from the region where qr — k is colhnear with £„. On the other hand Rl and 

come from the region where k is colhnear with £+. We see that the contributions Rl and Rl 
cancel with the corresponding contributions of diagrams 7a and 7b respectively. 

Contributions R^^ cancel with analogous contributions from diagrams 9al, a2 and 9bl, b2 
respectively, which have the form 

A ^^^(gal coUinear) 1 j^i l'cn\ 

c 



^ _^|!^(,„£t_i)„.-H|fli (60) 



1 



* ^^^(gbl coUincar) " D* D« fai\ 

Ai gg-^* = - R_ - R2 (.bij 



where we have not used the 'colhnear' superscript for those diagrams which do not have soft 
singularities. These contributions contain further artificial non-factorizable singularities 
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They in turn cancel with the contributions of diagrams 9a3, b3. 



Taking the sum of all these contributions and adding to it the quark self-energy contribution 
(|45|) , which we may call collinear as well, we obtain 



Ti(collinear) 




As mentioned earlier, this answer could have been written 'by hand', since the expression in 
brackets is easily recognized as arising from the quark splitting function 





dzPgg(z) = \ I + I ] dZ^il + d-Z'^ 



We see that the verification of (|65D in the Feynman gauge is complicated by the presence of 
spurious singular terms which mutually cancel. Nevertheless we prefer to use this gauge so as 
to minimize the possible error in the subsequent evaluation of the non-singular terms. 

3.3 Total singular virtual contribution 

The total singular part of the virtual correction to the Reggeon vertex for 7* qq is given 
by the sum of the 'soft -|- central' contribution of (^) and the 'collinear' contribution of (pSf). 
It is 



p(0)i 



(66) 



4. Real gluon emission 

Now we consider the contribution to the photon impact factor from the quark-antiquark- 
gluon intermediate state. First, we have to find the 7*i? — > qqg effective vertex. In principle, 
we could use the same approach as for the 7*/? qq vertex. Namely, we could calculate the 
amplitude for the conversion of a virtual photon into a qqg state in the high energy collision 
of the photon with a quark, take the part of this amplitude with a colour-octet state in the 
t channel and determine the 7*i? —>■ qqg vertex by comparison of this part with the Regge 
form (^. Since we need the vertex only in the Born approximation, it is sufficient to calculate 
the conversion amplitude in the same approximation and to take the Regge form in the lowest 
approximation in the coupling constant. 

The expression for the vertex can be obtained immediately, if we represent it by the diagrams 
of Fig. 10 with the vertex factor igf^2/s for the interaction of the Reggeon with quarks, and by 
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the vertex factor 7^^^^ (^) for the interaction of the Reggeon with gluons. It is straightforward 
to check that this direct calculation gives exactly the same 7*i? qqg vertex as the method 
described above. 

In the following we denote the momentum, polarisation vector and colour index of the 
emitted gluon by k, e and c respectively. The Sudakov decomposition is 



fc2 

k = (3pi + —p2 + k±, 
sp 

The vertex 7*/? — > qqg has the following properties. 



k — k I 



(67) 



1. In the 'soft' region 



/? < 1, fc^ < /i^ 



(68) 



the vertex is given by the formula for accompanying soft bremsstrahlung, i.e. 



P2 

k.p2. 



(g|tT|g) 



V2 



{q\ff\q) 



2. In the 'central' region 

/3 < 1, fc^ > Pij^, 
the vertex factorizes into a product of two Reggeon vertices, i.e 



(69) 
(70) 

(71) 



where g^' = qR + k, and where 'Jjiiqji', qn), the effective vertex for gluon production in 
a Reggeon- Reggeon collision M, is given by 



-qR'^ - qR^ + {k^ - 2ql,)^ - {k' - 2g^)| 



(72) 



Note that in ( |69D the gluon momentum is neglected in the vertex Tqq^^* . The neglect is valid in 
the whole region (|68[), not only when 2(k.£±) <^ /i^, (i.e. when k'^ <^ 13^^) where it is obviously 
true. In fact, it is a consequence of Gribov's theorem [|T^ concerning the region of applicability 
of formulas for accompanying soft bremsstrahlung. Thus the factorization (|69|) in the region 
(|68|) has a general nature. The same can be said about (|7lD . This equation was also checked by 
direct calculation; but, in fact, it also has quite a general nature and is connected with gluon 
Reggeization in QCD. 

The important point is that the soft ( |B5D and central ( 1701) regions overlap and cover the 
entire region (3 < xq with xq <^1. Let us denote 



1 



^ r,,-,,^.(g^)r^_^_^,,(g^) 



1191 



(73) 
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c.f. (18). The contribution of the qq g intermediate state to the impact factor is given by the 
integral of the function fqqg over 



W9 



X 



dx+ dx. d(3 d^-H+ d^-H- d^-'^k 



2x+ 2x_ 2(3 (27r)^-i (27r)^-i (27r)^-i^ 
c.f. ([T9|). In the soft region we find, using (B^), that 



(74) 



fqqg ~ fqq{^±^ ^±'i Qr^ Qr) ^ 



X 2g^N, 



c 



x^. 



x_ 



p%x+£^ - x_£+) 



{x+k - (31^ 



(75) 



|2 ^ (x_fc - [3£.f Nl {x+k - f3£+y{x.k - pL 
On the other hand in the central region it follows from (|7l|), ([7^) 

D-l 



where fqq is defined by (|T 
that 

2(27r)^"^ 

fqqg = fqq{x±, £±, Qr', Qr^) ^"7^ K ^°\q j^, , Q j^] q) 



Qr'Qr' 



(76) 



where q^, = qR + k; q'j^, = qj^, — q and K^^^ is the part of the non-forward LLA kernel 
entering (^2|) and defined by (p3|). Comparing ( [75| ) and ([76| ) one discovers an interesting result. 
The region of the applicability of the central form ( [7B[ ) for fg^g is larger than (^) and, in fact, 
is given by 

/? < 1, fc' > (3^fi\ (77) 

It is convenient to consider the contribution of the central region together with the subtrac- 
tion term in (|22|) . This term has two types of infrared singularities: the singularities contained 
in the trajectory function uj{t), and the singularities coming from the integration over q^^,. Ow- 
ing to properties (pT]) of the impact factor, the latter singularities are only due to the singular 
behaviour of the real part of K^^\qj^,, qr^; q) at qj^, — qr^ = = 0. It is useful to note that 
in (p2D In(so) appears as a factor of the total non-forward BFKL kernel. 



K{qR>,qR]q) = qRqR 5{qR, - q^) (^(t/j) + uj{t'^)) + Kr{qR, , qR, q), 



(78) 



taken in the Born approximation. Since the total kernel is infrared safe, the infrared singularities 
do not depend on the value of sq, so that we can take any suitable choice of sq to analyse them. 
It is convenient to take sq sufficiently large so that sq ^ /i^, but such that at the same time it 
satisfies sq <^ fJ^'^/xl- 

Then we can evaluate the integral over fggg in the region /? < xq as the sum of two subregions, 
in the first (second) of which k^ is less (larger) than (3'^so. The contribution of the first region 
can be calculated with the help of ([75|) . Its singular part is 



/ 



d$^°\^.(q«,q) 2f 



1 Iw 2 N 

+ - ln(xoSoJ 



1 1, 

^ + - In 



'xliL 



(79) 
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where 



(0) 

7* '.7' 



dM, 



271 



dp, 



qqi 



with the function fqq defined in (|T8|) and the phase space element in (|1^), so that 



d^ 



(0) 

7* '.7' 



7 ,7 



(80) 



(81) 



see (0) and ([T8|) . The contribution of the second region can be calculated with the help of 



(76). As already mentioned, it is convenient to consider this contribution together with the 
subtraction term in (p^). Our choice of Sq makes such a consideration especially convenient 
because of the total cancellation of the singular part of this contribution with the part of the 
subtraction term involving K^J^\ So, the singular contribution from the second region and the 
subtraction term is simply 

«?'..,•(<;.. 9)^ 'nf%^V (82) 



The dependence on sq vanishes in the sum of (|79| ) and (^2|), as it must. 



It remains to consider the singular contributions from the region (3 > Xq. These contribu- 
tions have a coUinear nature. For real emission we can use, without uncertainty, a physical 
gauge. In such a gauge, the 7*-R — > qqg vertex in the region of quasi-collinearity of the gluon 
and the quark (antiquark) momenta is given by the diagrams with a gluon line attached to the 
external quark (antiquark) line. Fig. lOal, bl (Fig. 10a3, b3). They can be easily calculated 
by the quasi- real electron method In this way the singular contribution from the region 
/? > Xq is found to be 



f2C, 
e N, 



c 




z 



d'^\>,^,{qn,q)-^ 



2 In I I - 3 



Xq 



(83) 



This completes the study of the real gluon emission contribution to the photon impact fac- 



tor. The singular parts of this contribution together with the subtraction term in eq. ( ^2]) are 
given by (^, and (jg^). 



5. Cancellation of the infrared singularities 

Let us denote the singular part of the virtual photon impact factor as the sum 



/ ,y* (^^, Q)singular 



rf<l>;°\,.(q^,g) 



rsoft + central _i_ rcoUinear _i_ rsoft _i_ ^central _i_ rcoUinear 
"virtual "virtual "real "real "real 



(84) 
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Here the first two terms in the square brackets are the singular parts connected with the virtual 
corrections to the 7*i? gg vertex and coming from the central and soft regions ( P9| ) and the 
coUinear regions 



rsoft + central 
'^virtual 



i5) 



rcoUinear 
'^virtual 



r 2c, 



e N, 



c 



2 In 



X-lX- 



Xn 



(86) 



The next three terms are the singular parts of the contribution of the real gluon emission taken 
together with the subtraction term in (^2]). They are given by (^), ( P^ ) and 



rsoft 
-^real 



r central 
■^real 



2r 



In 



- + - ln(a;oSo) 



^ + - In 



(87) 
(88) 



and 



6: 



coUinear 
real 



e Nr 



We see that the sum of the contributions (|8 
safe. 



2 In 



. ^0 , 



-3 



(89) 



(|89|) is zero. That is the impact factor is infrared 



6. Infrared safety of impact factors in the general case 

In the preceding Sections we have explicitly calculated the infrared contributions to the 
impact factor of the virtual photon coming from various sources, and have demonstrated their 
cancellation. In this Section we present general arguments, applicable to any impact factor, 
to prove its infrared safety. These arguments are based on well known theorems |^ about 
the absence of infrared divergences in totally inclusive quantities. Though these theorems are 
commonly used for cross sections and for the imaginary parts of forward scattering amplitudes, 
they apply equally well to non-forward scattering. 

Therefore we could be sure of the infrared safety of the impact factors provided that 
Reggeons were a kind of off-mass-shell partons, described by some fields, and the impact factors 
were the imaginary parts of the amplitudes for the scattering of these 'partons' off particles. 
But this is not the case for several reasons. First, the impact factors contain the subtraction 
(see (p2D where we have singularities due to the gluon trajectory function uj{t) as well as arising 
from K^^^f). Second, the region of integration over the phase space volume in ( [2^ ) is restricted 
by a theta-function. And finally, the Reggeon-particle vertices F* in (^2|) are not true field 
theory amplitudes. 

The first two reasons are evident and do not need further explanation. Let us amplify the 
third reason. We noted that at the Born level the vertices 7*i? qq and 7*i? — * qqg are given 
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by Feynman diagrams. It means, that these vertices appear as amphtudes of the effective field 



theory |T6|. An analogous statement is applicable in the general case, but it is only valid in the 
Born approximation. Recall that the corrections to the 7*-R qq vertex consist of two parts. 
One contribution, which comes from the central region of virtual gluon momenta (PP]), is given 
by ( ^9|) and is universal, i.e. process independent. The other contribution, which comes from 
the y4-region, is given by Feynman diagrams of the effective field theory with the restriction 
1/9 1 > on the Sudakov parameters of the gluon momenta. This last statement again has a 
general natur e, i.e. in the general case the one-loop corrections to Reggeon-particle vertices are 
given by the sum of the universal contribution (^9|) and a part given by the effective field theory 
with the restriction \P\ > /?o on gluon momenta. We can apply arguments about the absence 
of infrared divergences to the imaginary parts of Reggeon-particle scattering amplitudes in this 
theory. From the above discussion it is clear how impact factors differ from such imaginary 
parts. They differ by 



1) the existence of the subtraction term; 

2) the presence of the contribution of (virtual and real) gluons with \P\ < (3o; 

3) the limitation of the integration region over the phase space of the produced particles. 



So, to be convinced that the impact factors are indeed infrared safe, we need to check the 
cancellation of the infrared contributions to the impact factor coming from these three sources. 

It is convenient to use a notation analogous to (Q) for the singular part of the next-to- 
leading correction to the impact factor 



^A'A^qR, 9)singula 



A. _L ^virtual _|_ ^real _|_ 



(90) 



Then, taking into account the properties (|2T|) of the impact factors for colourless objects, we 
obtain from (|22|)-(g5|) 



5i 



(27r; 



D~l 



In 



+ u{tn) ln(:^j 



1 



-^M - + - In 



e2 ' e"^\q'lql^ 



Recall from Section 3 that the contribution of virtual gluons with 



the correction (^41) to Born vertex. Therefore, we obtain 



singular 

(91) 

< Po is universal. It gives 

(92) 



It is appropriate to consider the contribution of real gluons with < Po together with the 
limitation < s\ of the integration region over the phase space of the produced partons in 



20 



(22). This limitation does not play any role in the LLA, and becomes important only for gluon 
emission. Thus, in fact it is a limitation on the phase space of the emitted gluon to the region 



Recall that we take limit jSo before e ^ 0, so that for the calculation of the contribution of 



real gluons with \f3\ < (3o we can use the multi-Regge factorisation formula analogous to (|76D . 
Therefore we obtain 



±g'Nc 
(2' 



TT 



d'^-^k fPo dp 



- singular 



\ + - ln(sA/?o; 



(93) 



We see that the sum of (pl|)-(p^) is zero, which means that we have verified the infrared safety 
of impact factors in general. 



7. Discussion 

In this paper we have investigated, to the next-to-leading order accuracy, the infrared prop- 
erties of impact factors appearing in the BFKL description of small x processes. We have shown, 
by explicit calculation, the cancellation of the infrared singularities in the impact factor of the 
virtual photon. The virtual photon impact factor was considered because of the importance 
of corrections to LLA for diffractive qq electroproduction and because of the possibility of cal- 
culating this impact factor in perturbation theory for large enough photon virtuality. We also 
presented a general proof of the infrared safety of the impact factors describing the transitions 
between colourless particles. 

We have considered impact factors in the general case of non-forward scattering and have 
used the definition of impact factors given in [|], |^. These impact factors are defined for the 
case when the energy scale sq (in the Mellin transform factor {s/sq)^ in (^) is a constant 
independent of the integration variables (the Reggeon "masses"). One may prefer to use a 
scale which is dependent on these variables. (For instance in the particular case of forward 



scattering, i.e. for calculation of cross sections, the scale y^iqi was used in [Q). It was shown 
in [1^ that to the NLA accuracy one can indeed change the scale sq in (^ for any factorizable 
scale /1/2, with /j depending on (and on q), without changing the Green function, provided 
that the impact factors are also changed according to 

^A'A{qR,q]SQ) ^A'AiqR,q;so) 

^ J qR'qR' \soj 
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where K^^^ is the total non-forward LLA kernel, 



K^''\qR>,qR;q) = Ki^\qn> , Qr, q) + q'Ww 6{q^- q^,) {uj{tR) + ^(4)). (95) 

K^^^ is the part of the LLA kernel related to real gluon production, (|23|), and uj{t) is the gluon 
Regge trajectory. The properties of the impact factors (EO) and the infrared safety of the total 
kernel (due to cancellation of the infrared singularities connected with K^^"* and with the gluon 
trajectory) guarantees that the transformation ( plf ) does not change the infrared properties of 
the impact factors. 

We emphasize that the above definitions differ from those introduced, for the case of the 
forward scattering, in [|T3| and used in and in [Q, where impact factors of partons (quarks 
and gluons) were calculated for a colour singlet state in the t channel. The definition advocated 
in these papers is based on "the requirement that the subtracted leading term should satisfy 
at finite energies for the impact factors to be well defined (e.g., without spurious infrared di- 



vergencies)" |T^. In order to fulfil this requirement additional operator factors Hl^Hr were 
introduced 0, in the Green function G^. In our opinion, the only physical requirement 



is the absence of infrared singularities in the impact factors of colourless objects, and we have 
explicitly demonstrated here that our definition of the impact factors fulfills this important 
requirement. 
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Figure Captions 

Fig. 1 Schematic representation of the process A + B —>■ A' + B' displaying the structure of (H). 
The zig-zag hues represent Reggeized gluon exchange. 

Fig. 2 The intermediate process A + B ^ a + b mediated by Reggeized gluon exchange. 

Fig. 3 The lowest-order Feynman diagrams for the process 7*g — > {qq)q, showing the particle 
four momenta. 

Fig. 4 Schematic representation of the vertex r^^^J^,, showing, at lowest order, the transition 
7*i? ^ qq. 

Fig. 5 The two-gluon exchange Feynman diagrams mediating the process 7*g —>■ {qq)q- 

Fig. 6 Feynman diagrams mediating 7*g —>■ {qq)q which must be considered with those of Fig. 5. 

Fig. 7 Diagrams containing the gluon-gluon-Reggeon vertex which contain the corrections to the 
^qq,^* coming from Figs. 5 and 6. 

Fig. 8 The gluon-gluon-Reggeon vertex of (^T]). 

Fig. 9 Loop diagrams within the produced qq pair. The quark (antiquark) self-energy diagrams 
are not shown. 

Fig. 10 The real gluon emission diagrams from the ^qq,-y* vertex. 
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Fig.7 
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